DAHLGREN  DIVISION 

NAVAL  SURFACE  WARFARE  CENTER 

Dahlgren,  Virginia  22448-5100 


NSWCDD/TR-94/15 


AD  “  A278  876 


PRELIMINARY  STUDY  OF  SATELLITE  ATTITUDE 
DETERMINATION  FOR  THE  LANDSAT  7  SPACECRAFT 


BY  PATRICK  J.  FELL 

STRATEGIC  AND  SPACE  SYSTEMS  DEPARTMENT 


APRIL  1994 


Approved  for  public  release;  distribution  is  unlimited. 


a^ALITf  lUSPEOrBD  i 

94  5  06 


DTIC 

iELECTEf% 


NSWCDD/TR-94/15 


FOREWORD 


This  report  explores  the  use  of  multiple  data  sources  to  support  satellite 
attitude  determination  for  the  LANDSAT  7  spacecraft.  This  study  was  initiated  at 
the  direction  of  the  Defense  Mapping  Agency  (DMA)  as  a  means  of  understanding 
how  attitude  could  be  refined  in  order  to  exploit  LANDSAT  7  stereo  pairs  for 
mapping  and  charting. 

The  work  was  performed  in  the  Space  and  Surface  Systems  Division  of  the 
Strategic  and  Space  Systems  Department.  This  report  has  been  reviewed  by  James 
L.  Sloop,  Head,  Space  and  Surface  Systems  Division. 


Approved  by: 


3Y  L.  SCHMIDT,  Head 
Strategic  and  Space  Systems  Department 
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ABSTRACT 


This  report  presents  the  results  of  an  error  analysis  for  satellite  attitude 
determination  for  the  LANDS  AT  7  spacecraft.  The  purpose  of  this  study  was  (1)  to 
evaluate  the  contribution  of  dual-antenna  Global  Positioning  System  (GPS) 
measurements  to  spacecraft  attitude  uncertainty  and  (2)  to  determine  if  spacecraft 
attitude  could  be  determined  to  an  accuracy  that  would  allow  direct  mapping  from 
LANDSAT  7  stereo  pair  imagery.  Although  preliminary,  the  results  indicate  that 
GPS  measurements  have  the  potential  to  improve  spacecraft  attitude  if  certain 
technical  goals  can  be  achieved  associated  with  GPS  instrumentation.  The  useability 
of  stand-alone  stereo  pairs  for  direct  mapping,  satisfying  Defense  Mapping  Agency 
(DMA)  product  standards,  is  questionable  due  to  correlated  stellar  sensor  errors. 
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INTRODUCTION 


The  ability  to  estimate  precise  attitude  for  the  LANDSAT  7  spacecraft  will  be  a 
determining  factor  for  exploiting  single-pass  stereo  pair  imagery  for  mapping 
without  the  need  for  large  area  photogrammetric  triangulation.  As  initially 
designed,  spacecraft  attitude  during  the  mission  would  be  determined  from  a 
combination  of  stellar  sensor  and  rate  gyro  data.  However,  these  sources  of  orienta¬ 
tion  information  are  not  sufficiently  accurate  to  preclude  the  need  for  photo¬ 
grammetric  triangulation.  Therefore,  it  was  prudent  to  explore  alternative  means 
for  possible  enhancement  of  mission  attitude  determination.  The  data  sources 
considered  to  augment  the  design  were  Global  Positioning  System  (GPS)  dual¬ 
antenna  phase  measurements  and  stereo  pair  mensuration. 

As  a  first  step  in  understanding  what  accuracy  might  be  achievable,  a 
preliminary,  one-dimensional  spacecraft  attitude  study  was  performed.  This  study 
examined  four  data  types  in  a  least  squares  covariance  analysis  of  satellite  attitude 
determination.  The  measurement  types  considered  are  listed  in  Table  1 .  The  rate  at 
which  these  observations  would  be  acquired  or  processed  and  their  nominal  precision, 
in  terms  of  how  accurately  orientation  can  be  obtained  from  the  measurements,  is 
also  provided.  Nominal  precision  should  be  interpreted  to  mean  subsystem 
specifications  in  the  case  of  stellar  attitude  and  rate  gyro  stability,  a  technology 
performance  goal  for  GPS  attitude  determination,  and  an  anticipated  capability  from 
stereo  mensuration  based  on  image  pixel  size  and  ground  coverage  provided  by  a 
typical  stereo  pair. 

TABLE  1.  LANDSAT  7  DATA  TYPES  USED  IN  COVARIANCE  ANALYSIS 


Data  Type 

Data  Rate 

Nominal  Precision 

Stellar  Attitude  (Celestial 

Sensor  Assembly) 

2  per  image  operation 

4  arcsec  random  bias 

2.7  arcsec  random  noise 

Gyro  Derived  Attitude  Change 

1  per  second* 

0.01333  arcsec/sec 

GPS  Attitude  (Phase  Derived) 

1  per  second 

36  arcsec2 

Stereo  Pair  Mensuration 
(Relative  Orientation) 

1  per  image  pair 

5.8  arcsec 

10. 1  sec  rate  is  available,  1  sec  rate  considered  in  analysis 

2Based  on  goals  of  a  Naval  Research  Laboratory  (NRL)  Satellite  Program  (1995) 
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Figure  1  provides  a  sense  of  the  geometry  of  LANDSAT  7  stereo  pair 
acquisition.  The  imagery  and  support  data  collection  takes  place  during  a  3-min 
interval.  This  image  acquisition  scenario  was  the  basis  for  a  covariance  analysis, 
which  assumed  the  data  sources  described  in  Table  1  were  collected  or  derived  from 
measurements  at  the  specified  rates.  Stellar  attitude  was  assumed  to  be  measured  at 
the  endpoints  of  the  time  interval  spanning  Ti  to  T4.  The  images  were  assumed  to  be 
instantaneously  acquired  at  T2  and  T3  separated  in  time  by  120  sec.  In  addition  to 
stellar  attitude,  orientation  changes  were  assumed  to  be  continuously  provided  by 
rate  gyro  data  and  it  was  assumed  that  GPS  orientation  was  determined  from  dual¬ 
antenna  phase  measurements  at  a  1-sec  rate.  Image  mensuration  provided  the 
relative  orientation  of  the  stereo  pair. 


FIGURE  1.  GEOMETRY  OF  LANDSAT  7  IMAGERY 
(OVERLAP  UNDERREPRESENTED) 
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ATTITUDE  DYNAMICS  MODEL 


Over  the  time  interval  from  Ti  to  T4,  during  the  image  operation,  the  attitude 
dynamics  of  the  spacecraft  (one-dimensional  model)  was  assumed  to  be  given  by  a 
polynomial  of  degree  three: 


0  =  at^  +  bt^ ct  -l-  d 

where  a,  b,  c,  and  d  are  model  parameters  to  be  estimated  and  t  is  the  number  of 
seconds  from  Ti.  Of  particular  importance  are  the  uncertainties  in  spacecraft  atti¬ 
tude  at  imaging  times  T2  and  T3. 

Since  stellar  and  GPS  attitude  provide  absolute  measures  of  orientation,  these 
data  will  contribute  to  all  attitude  model  unknowns.  Rate  gyro  data  and  image  men¬ 
suration  will  provide  angular  change  and  relative  image  orientation,  respectively, 
and  thus,  contribute  to  a  determination  of  only  three  of  the  four  parameters  of  the 
attitude  model. 


MEASUREMENT  MODELS  AND  STATISTICAL  PROPERTIES 


This  section  presents  the  measurement  models  that  were  incorporated  into  the 
spacecraft  attitude  analysis  presented  below.  For  each  measurement  type,  the 
mathematical  model  is  provided  along  with  the  measurement’s  statistical  repre¬ 
sentation  and  the  design  matrix  used  in  the  development  of  the  covariance  analysis. 


STELLAR  ATTITUDE 

Stellar  attitude  data  were  assumed  to  be  measured  at  each  end  of  a  180-sec 
imaging  interval  assuming  an  error  structure  consisting  of  a  random  bias  and  a 
random  (white  noise)  component.  At  times  Ti  and  T4,  the  stellar  sensor  will  provide 
estimates  of  spacecraft  orientation  given  by 

e‘’^  =  0^  +  e(Tj) 

and 

0°  =  0  +e(T  ) 

4  4  4 
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where  ei  =  e  (Ti)  and  04  =  e  (T4)  represent  errors  in  stellar  attitude  having  the 
characteristics  (see  Table  1): 

=  a  +  Hi 

£4  =  a  + 

where  the  random  bias  a  has  a  nominal  standard  deviation  of  4.0  arcsec  and  the 
random  error  q  has  a  nominal  standard  deviation  of  2.7  arcsec.  Thus,  the  variance 
and  covariance  for  the  errors  ei  and  02  are  given  by 

E  =E  =a‘^  +  E  q^q^  =a^  +  E  q^q^ 

The  covariance  matrix  for  the  pair  of  stellar  updates  is  thus, 

+  E  q^q^ 

a'^  +  E  q^q^ 

which  yields  an  error  correlation  of  approximately  0.7  (correlation  coefficient).  The 
inverse  of  this  matrix  provides  the  weighting  for  these  measurements  when  used  in  a 
least  squares  covariance  analysis. 

For  stellar  attitudes  at  times  Ti  (0.0  sec)  and  T4,  the  design  matrix  for  the 
covariance  analysis  is  given  by 

0  0  0  1 

A  = 

STELLAR  rp  fp  2  fp  « 

.  4  4  4  ^ 

GLOBAL  POSITIONING  SYSTEM 

The  GPS  measurements  (phase)  are  assumed  to  be  processed  to  obtain  space¬ 
craft  orientation  every  second  within  the  interval  (Ti  T4]: 

0“  =  e  +  V 
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where  the  errors  v  were  assumed  to  be  statistically  independent.  Accuracy  levels  of 
36  arcsec  (NRL  program  goal  for  1995),  180  arcsec  (demonstrated  accuracy  at  the 
Naval  Surface  Warfare  Center  Dahlgren  Division  (NSWCDD))  and  500  arcsec  were 
considered  in  the  covariance  analysis  discussed  below.  Thus,  the  covariance  matrix 
for  the  set  of  GPS  derived  orientations  is  given  by 


flvQVol  0  0 . 0 


S  = 


0  0 . 0 

0  0  0  £lv  V  1 

n  n 


The  design  matrix  for  the  GPS  derived  orientations  is  given  by 


^GPS 


0  0  0  1 


t  t  ^  t  1 
‘2  *2  *2  ^ 


I  f  ^  f  ^  t  i 

L  n  a  n  J 

where  is  the  number  of  seconds  from  Tj  (taken  to  be  the  time  reference  epoch  in  this 
analysis). 


RATE  GYRO 

The  rate  gyro  provides  (by  integration)  the  change  in  0  over  any  time  interval 
[tj,  t2Y 

f '2“  ( '2'  <-('2  -  'll  +  V,, 

Although  available  at  a  higher  rate,  it  was  assumed  that  the  gyro  data  would  be 
sampled  at  a  1-sec  rate  from  Ti. 

The  error  y  in  A0  is  the  integrated  effect  of  the  gyro  drift  rate  and  thus, 
represents  a  correlated  error  signal  throughout  any  3-min  interval.  Two  statistical 
gyro  error  models  were  considered.  The  first  model  represented  constant  drift  noise 
with  a  0.01333  arcsec/sec  uncertainty,  while  the  second  model  consisted  of  drift  noise 
and  drift  noise  rate: 
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o.  =  0.01333  (Modell) 
t-o 

o  .  = 

0) 

u  =  0.01333  {Model2) 

V  =  0.000079 

where  the  units  are  arcsec/sec  and  t  is  in  seconds  from  T;  when  the  first  stellar 
update  is  made.  These  error  models  were  taken  from  LANDSAT  7  technical 
materials  provided  by  DMA  and  represent  the  uncertainty  in  gyro  drift  at  the  one  (1) 
sigma  level. 

Given  that  is  the  integrated  output  of  the  rate  gyro  over  [ti,  t2l,  then  the 
uncertainty  in  the  error  Yi2  (assuming  ti  =0=Ti)  is 


f '2 

=  a.  dt 
Jo  eo 

=  0.01333<.^  (Model  W 

=  (2/3)u<2^+(2/5)v<2^^  (Model2) 


For  two  distinct  times  t2  and  ta,  the  errors  in  AB\2  and  will  be  correlated 
depending  on  the  separation  of  t2  from  ta.  The  closer  t2  is  to  ta  the  higher  will  be  the 
correlation  since  the  errors  as  manifested  by  gyro  drift  become  more  identical. 

For  t2  and  ta 


Yu  =  Yi2  +  dt 


1::  • 

Multiplying  through  by  Y12  gives 


Yu  Yt2  Yi2  ¥^2  Yj2 


(*3  . 

^dt 


Taking  expectations  yields 


^Yi3Yi2l  =  ^Yi2Yi2l 


'12 


^dt 
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The  second  term  on  the  right-hand  side  of  this  last  equation  is  the  covariance 
between  the  accumulated  error  in  the  gyro  alignment  Y12  up  to  t2  and  the  additional 
error,  which  accumulates  between  t2  and  ts.  For  the  purposes  of  this  analysis,  it  was 
assumed  that  these  errors  would  be  uncorrelated.  Thus,  the  covariance  between  the 
errors  Y12  and  Y13  was  taken  as  the  variance  of  Y12.  This  assumption  essentially 
implies  that  the  errors  in  gyro  drift  can  be  treated  as  a  random  walk.  This 
assumption  may  be  weak;  however,  no  detailed  performance  data  on  the  gyro  types 
being  considered  for  LANDSAT  7  were  available  to  refine  this  assumption. 
Accordingly,  the  variance-covariance  matrix  for  orientation  changes  after  Ti  based 
on  sampled  rate  gyro  data  is  given  by 


Sy 


^Yi2Yi2'  ^Y,.,Y,.,J 

^Yi,Yi,I  £lYi3Y,3l 


£Iy„y„1 

^Y,3Yi3l 


^Yi2Yi2l 


®Y,3Y,3l 


^Y„„Yi,„I 


The  design  matrix  for  gyro  derived  orientation  changes  AO  from  ti  =  0  =  Ti  to 
times  ti  (i  =  2, . m)  is  given  by 


A 

GYRO 


t.  0 

o 


The  times  ti  occur  each  second  after  the  first  stellar  update  at  time  Ti  and  end  with 
the  second  stellar  update  at  T4. 


IMAGE  MENSURATION 

It  was  assumed  that  image  mensuration  at  points  separated  by  250  km  would 
provide  relative  orientation  at  image  times  T2  and  T3  to  a  nominal  accuracy  of 
5.8  arcsec.  This  uncertainty  was  based  on  a  mensuration  error  of  1  pixel  (5  m).  The 
observation  model  for  image  pair  mensuration  is  given  by 

A0  =  0.,— 6.,  +  a 
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where  02  and  83  are  the  spacecraft  orientations  at  image  times  T2  and  T3  and  a  is  the 
error  in  relative  orientation  derived  from  the  mensuration  process.  The  variance  of 
the  relative  orientation  A0  is  given  by  E(aal. 

The  design  matrix  based  on  this  observation  is  given  by 


^MENSURATION 


*^2  _  'P  ^ 

3  2 


0 


COVARIANCE  ANALYSIS 


Using  these  data  structures,  data  rates,  and  noise  processes,  least  squares 
normal  matrices  were  developed  for  each  of  the  four  observation  types  consistent  with 
the  image  acquisition  scenario  of  Figure  1; 


where  for  each  observation  type  the  weight  matrix  Wi  is  the  inverse  of  the 
corresponding  variance-covariance  matrix  for  that  observation  type  as  given  above. 
After  formation,  the  matrix  ATWA  was  inverted  to  produce  the  variance-covariance 
matrix  for  the  determination  of  the  attitude  modeling  parameters  a,  b,  c,  and  d. 
Using  this  result,  the  uncertainty  in  spacecraft  orientation  at  any  time  during  the 
interval  can  be  developed  through  linear  error  propagation  using  the  following 
equation: 


E 


0(n 


^  e  t  I 


A^WA 


1 


T 


In  particular,  the  uncertainty  in  6(t)  at  image  times  T2  and  T3  were  of  interest. 

The  results  from  this  covariance  analysis  are  given  in  Table  2  for  various  cases, 
which  considered  nominal  measurement  precision  as  well  as  certain  variations  in 
order  to  gain  insight  into  the  influence  of  each  error  source.  Case  1  was  based  on  the 
use  of  the  nominal  error  statistics  from  Table  1.  Cases  2  through  7  provided  orienta¬ 
tion  results  as  a  function  of  variations  in  the  noise  processes  associated  with  the 
measurements. 


Several  conclusions  can  be  drawn  from  the  results  of  this  preliminary  analysis: 

(1)  Simple  stereo  pair  mensuration  has  little  impact  on  spacecraft  attitude 
accuracy  (Case  2  vs  Case  4)  since  the  correlated  error  in  orientation  of  a  stereo  pair 
due  to  the  stellar  sensor  cannot  be  reduced  through  mensuration. 
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(2)  The  variation  in  the  gyro  noise  models  considered  here  had  no  significant 
impact  on  the  overall  results  (see  Cases  1, 2, 6,  and  7). 

(3)  Attitude  accuracy  during  the  imaging  operation  appears  to  be  a  linear 
function  of  stellar  sensor  accuracy  (see  Cases  2  and  5). 

(4)  GPS  attitude  determination  worse  than  180  arcsec  has  little  impact  on  the 
results  (Case  2  vs  Case  3).  If  GPS  attitude  is  improved  significantly,  to  a  level 
consistent  with  an  NRL  technology  goal  of  36  arcsec,  then  spacecraft  attitude  at  the 
imaging  times  could  be  significantly  improved  (Case  1  vs  Case  2  and  Case  7  vs 
Case  6). 


TABLE  2.  COVARIANCE  ANALYSIS  RESULTS  FOR  LANDSAT  7  IMAGE  PAIR 

(180-SEC  IMAGE  OPERATION) 


Case 

GPS 

Stellar 

Gyro 

Mensuration 

Attitude  Result 

1 

36 

4.0/2, 7 

0.01333 

5.8 

2.3 

2 

180 

4.0/2.7 

0.01333 

5.8 

4.2 

3 

500 

4.0/2.7 

0.01333 

5.8 

4.4 

4 

180 

4.0/2.7 

0.01333 

2.9 

4.2 

5 

180 

2.0/1.35 

0.01333 

5.8 

2.2 

6 

180 

4.0/2.7 

0.01333/1**0.5 

-1-0.000079/1**1.5 

5.8 

4.2 

7 

36 

4.0/2.7 

0.01333/t**0.5 

-l-0.000079/t**1.5 

5.8 

2.3 

Note:  Attitude  results  are  in  units  of  arcseconds 


SUMMARY 


Although  limited  in  scope,  this  analysis  has  provided  an  indication  of  the 
impact  that  certain  measurement  types  would  have  on  determining  the  attitude  of 
the  LANDSAT  7  spacecraft.  It  is  clear  that  high  precision  attitude  will  depend 
critically  on  the  accuracy  of  the  stellar  sensor  system.  Augmentation  with  GPS  may 
prove  useful  if  certain  GPS  technical  goals  can  be  achieved.  However,  the  ability  to 
exploit  (in  a  metric  sense)  LANDSAT  7  imagery  for  mapping  and  charting  using 
stand-alone  stereo  pairs  may  be  marginal  considering  product  standards  for 
positional  accuracy  and  the  results  from  this  analysis. 
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